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INTRODUCTION 

Many papers have been written on magnetic, bearings. These papers are 
divided into two categories: One is the use of magnetic bearings taking the 

place of contact bearings, and the other is the use of noncontact actuators 
for active vibration control. These catergories are combined in applications. 
The research on using magnetic bearings in place of contact bearings has been 
going on for years. As early as 1842 it was known that magnetic bearings 
should have at least one direction .actively cpntrol led .(.ref . 1). For the last 
century and a half there has been a great deal of research and development. 
Magnetic bearings have been used mainly in experiments and in aircraft and. mis 
sile inertial guidance control. However, these bearings are now coming into 
use in the industrial world (ref. 2) for turbo-molecular pumps (ref. 3), mach- 
ine tool spindles, etc. Although these magnetic bearings are for rigid rotors 
(refs. 4 to 6), they recently have been used with flexible rotors (refs. 7 to 
9). In such applications, the magnetic bearings suppress unbalance vibration. 
Many papers have been published on this subject (refs. 10 to 16). 

This paper is concerned with active vibration control of flexible rotors ' 
supported by magnetic bearings. Using a finite-element method for the mathe- 
matical model of a flexible rotor, the author has formulated the eigenvalue 
problem taking into account the interaction between the mechanical system of a 
flexible rotor and the electrical system of magnetic bearings and a controller 
However, for the sake of simplicity, gyroscopic effects are disregarded herein 
It is possible to adapt this formulation to a general flexible rotor - magneti 
bearing system. Controllability with and without collocation sensors and actu 
ators located at the same distance along the rotor axis is discussed for the 
higher order flexible modes of the test rig. In conclusion, the author has 
proposed that it is necessary to add new active control loops for the higher 
flexible modes even i n the case of collocation.' It is then possible to stabi- 
lize for the case of uncol location by means of this method. 


ANALYSIS OF ROTOR - MAGNETIC BEARINGS SYSTEM FOR RIGID ROTOR 

Basic Equation 

This paper deals with a radial-type magnetic bearing system. Figure 1 
shows four assembled electromagnets. The whole sys.tem is shown in figure 2 
where the flexible rotor shown is rigid. We have the following assumptions: 

(1) The attractive forces are proportional to the square of the coil 
current. 

(2) Both electromagnetic bearings have the same, characteri sties . 
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(3) , The. induced voltages of electromagnets are ignored. 

(4) The coil inductances are independent of frequency and gap length and 

are constant. ...... .... ' 

(5) This system is uncoupled between the x_ and y directions. 

( 6 ) Only small vibrations near equilibrium are considered. 

(7) The rotor mass is concentrated at the bearing. 


Under these conditions, we obtain the following equation for one electro- 
magnet:'; ‘ ' 


E - L Hi * RI 


( 1 ) 


where E is the coil input voltage, L the coil inductance, R the coil 
resistance, and I the coi 1 ; current. The attractive force of an electromagnet 
cain generally be given by 


P = 


p 0 an 2 i 2 

H 2 


( 2 ) 


where P is the attractive force, po the permeability, A the face area, N 
the number of, winding turns, and H the gap length. From the standpoint of 
small vibration near equilibrium, P, H, and I are given by 

p s P 0 ♦ P 


where po is 'the steady-state attractive force, io the steady-state current, 
ho the steady-state gap length, p the control . attractive force, i the con- 
trol current, and h the control gap length. Using the Taylor series expan- 
sion for small values of i and h and assuming io >>i and ho >>h, we can 
finally get the following attractive control force in linear terms: 





(4) 


The four control attractive forces Pi , P 2 , P 3 , and P 4 are shown in figure 1. 
The resultant forces of the horizontal control attractive force and the verti- 
cal control attractive force are, respectively, 


mx = p 1 - p 3 

my- = <P S .+ P 2 > - <P Q + P 4 ) - mg 


(5) 
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where m is the half-mass of the rigid rotor and g the acceleration of grav- 
ity. Assuming that control currents to coils are i 3 = -i] and i 4 = - i 2 , we 
have the following expressions considering equation (4): 



where the control gap length h is replaced by the displacements x and y 
from equilibrium; mg + pq is replaced by the steady-state attractive force 
p s and the steady-state current i s . 


The phase lead circuit shown in figure 3 is used to compensate for the 
time constant of electromagnets in order to maintain control at higher frequen- 
cies. Using figure 3, we can obtain the following equation: 


dV E 1 R 1 * R 2 
dt ‘ CR ] CR ] R 2 


(7) 


where the new parameter V is the voltage at R 2 and C. 


State Equation and Control System 

The integral values of displacements are used as control variables to com- 
pensate for any drift from the equilibrium position. The formulas are 



Replacing E with Ei - V in equation (1) (since E = E] - V), we get 


dj_ 1 
dt = L 


(E ] - V) 




Using equations (5) to (9), we can write the state equations as 

{X} = CA x ]{X} + {b}u x 

{X} = [Ay] {Y} + {b}Uy 


(9) 


( 10 ) 

( 11 ) 


where [A x ] is the horizontal system matrix and [Ay] is the vertical system 
matrix. There are as follows: 
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[A x ] = 


0 ' • 1 o • ... 0 


A h 

mh. 


0 0 


CAy] = 


10 0 

1 op 

0 0 0 

1 0 


+ 0 0 


; 4Pc 

mi. 


R 

L 


0 

0 


+ R 2 

cr 1 r 2 


( 12 ) 


m l h 


0 Oy 
1 


4 1 i . ' i. 
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0 0 
0 0 

0 0 


s "0 

0 

R 
L 


\ 


0 

0 

0 

1 

L 


R.| + R 2 
CR 1 R 2 J 


( 13 ) 


The state vectors {X} and {Y}, the control vector {b}, and the control 
inputs u x and Uy are given by 


II 

X 

w-' 

1 1 

X 

X • 

p 

X 

-< 

II 

y. y, 3 . iy, 

II 

JO 

— l-i 
O 
O 
O 

1 1 


u x “ E 1 x | 


C 

*< . 

II 

m . 
<< 


( 14 ) 

( 15 ) 

( 16 ) 
( 17 ) 


where T means transpose. The control input voltages E] x and Ely are 
given by combining the feedback wi th .displacements , velocities, and integra- 
tions of displacements. 
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Therefore, the output feedback control system of PID type is 


y - - (F y w> 

where 

< F x> = ,[ f lx’ f 2x’ f 3x’ °’ 

= [ f ly * f 2y’ f 3y’ °’ 

Elements of { F x } and {Fy} are feedback gains. 

Characteristic Equation 

Substituting equation (18) into equations (10) and (11) and taking the 
Laplace transformation give the characteristic equation In the horizontal 
direction as 

5 4 3 2 

s + a^s + a 2 s + a 3 s + a 4 s + a 5 = 0 (20) 

where • •- 

a, - £ * Q 

. OR 4 Pp f 2 % 

a 2 “ L migL " mhQ 

4 P 0 f 2° 4 P 0 f l 4 Pp R 4 PpO 4 Pp f 2 

a 3 mi -L ml^L mh^L mh- + mi.LCft^ 

4 Pp f 3 4 Pp f lQ 4 Pp R Q 4 Pp f l 

■ ■ a 4 mi-l ' ml qL mii-L * mi-'.CR. 

4c P f 3 C 4 Pp f 3 
a 5 " ‘ ml q L * mt 0 LCR, 


Likewise,. in the .verti cal direction, 

s 5 + b^ 4 + b 2 s^ + b 3 s 2 + b 4 s + b 5 = 0 (21) 




>] 


( 18 ) 


(19) 
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where 


b i -T + Q 

b, V'.AQ.a 


f , B f.BQ 

b 3 = "I” + ~T~ ~ L 


AR f 2 B 


LCR, 


- A Q 


a V + _ ARQ _ ^ 


LCR, 


f 3 BQ f 3 8 


LCR, 


A . 1 f !i + V 

m \^h Q h Q; 


m Vs ’0; 


Q = 


R^ + Rg 

cr 1 r 2 


The subscripts x and y of feedback gains are omitted after this section. 
The stability criterion of feedback gains f i , f 2 , and f 3 is decided by 
equations (20) and <21 ) ■ ' 


ANALYSIS OF ROTOR - MAGNETIC BEARING SYSTEM FOR FLEXIBLE ROTOR 

The mathematical , model for this flexible rotor is given by using a finite 
element method: • • • 


CM] {q} + CGl.{q} + CK] {q} = {0} (22) 

where CM] , CG], and CK] are the mass, gyroscopic, and stiffness matrices, 
respectively. We assume CG] * 0 for simplicity. This flexible rotor is now 
divided into six elements as shown in figure 2, and the state values are set 
at X] and e, at each station. CM] and CK] become 14x14 matrices and {X} 
vector has a length of 14 in this case. The flexible rotor is supported by 
magnetic bearings at stations 2 and 6, and sensors are located at stations 1 
and 7. From these results, the formulations in : the horizontal direction of 
the left magnetic bearing shown in figure 2 are defined by 
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4P 4P 

U x& " H a X 2 I a 1 xQ. 


E 1xH ' V x4 - L d af - R V 


V 


’1 X2, 


xft " CR 


R 1 + R2 
CR ] R 2 


XJL 


( 23 ) 


E i xa = f u x i ” f 2a x i f 32, a a 


where the force u x q, acts on station 2, and the subscript a indicates the 
left magnetic bearing. The steady-state values po, hQ, and io are replaced 
by Pq,, Hq,, and Iq,. The expressions are derived about the right magnetic 
bearing in the same manner. Combining equations (22) and (23), we can obtain 
the 34th order state equation as shown in equation (24): 


[See page 8] 
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-[M] - [K] 
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element (17,3) *6: element (25,11) *11: element (31,1) *16: element (29,19) *21: element (34,9) 

element (17,29) *7: element (25,32) *12: element (31,15) *17: element (31,5) *22: element (34,23) 

element (30,1) *8: element (33,13) *13: element (34,13) *18: element (31 ,19) 

element (29,1) *9: element (32,13) *14: element (34,27) *19: element (32,9) 

element (29,15) *10: element (32,27) *15: element (29,5) *20: element (32,23) 














where the subscript r means the right magnetic bearing. In addition, 
fj fj r , and f^ show feedback gains added for active vibration control 

loops for higher modes. Similarly, 'we get almost the same equation for the 
vertical direction. Equation (24) is systematically formulated for a general 
flexible rotor with magnetic bearings. We can design the control system by 
solving the eigenvalue problem of equation (24). 


• EXPERIMENTS 
Test Rig 

The laboratory test rotor is shown schematically in figure 2. Figure 4 is 
a photograph of the test fig. The flexible shaft, which is a 10-mm-diameter 
drill rod that is 1000 mm long, has five disks with masses as indicated in 
figure 2. The total mass of the flexible rotor is about 3.3 kg. Two of the 
disks are used for magnetic bearings and two others are used for sensors - i.e., 
the bearings and sensors are not collocated. The rotor is supported by two 
magnetic bearings with a bearing span length of 680 mm. The shaft is connected 
to the driving. motor by a flexible coupling. This is a radial magnetic bearing 
system with four degrees of freedom only because the radial directions are non- 
contacting. The gap between the electromagnets and the disk is about 2.0 mm. 

The backup beari ngs- protect in an emergency when the amplitudes exceed 1.0 mm. 
The shaft displacements, which are measured by gap sensors of the eddy current 
type, are measured as close as possible to the magnetic bearings. These sig- 
nals are the inputs for the analog controller. The analog controller consists 
of proportional , differential , and Integral circuits, the phase lead circuit, 
and some filters. The control signals from the analog controllers are supplied 
to the electromagnets through servoampl ifiers . 


. Experiments for Rigid Rotor 

At low frequencies the test rig rotor can be considered rigid. Figure 5 
shows the three-dimensional stability criterion chart for the horizontal direc- 
tion based on the Hurwitz's criterion of equation (20). Figure 6 is for the 
vertical direction. It is clear from the results that feedback gains f ] , f^, 
and f 3 should be negative. The stability improves when the absolute values 
of feedback gains become large; however, the stability region becomes narrow 
when | f i | is large. 

Figure 7 shows the root loci for different feedback gains. . The calcula- 
tions shown were done on the basis of feedback gains at the intersection 
points of the solid line, the dotted lines, and the dot-dash lines. In the 
case of the solid line (f i ) , the system becomes unstable because the pole which 
goes to the right and up passes over the imaginary axis to the other side. 

This corresponds to the case in which the stabi 1 ity region becomes' narrow when 
| f 1 1 becomes large (see figs. 5 and 6 ). In the case of the dashed line (f 2 >. 
the system does not become unstable though one pole goes toward the origin. 
Figure 7 shows that the sensitivity for variation of f .3 is very low when the 
behaviors of real roots are omitted^. The static spring constants were measured 
with the rotor floating and the system statically stable. This was done by 
measuring the deflection of the' shaft and the force of the tension spring bal- 
ance after attaching the spring 1 balance to the end of the shaft. The results 
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are shown in figure 8. The force vs. deflection curves for these magnetic 
bearings are linear over the range of interest. The spring constant used in 
the following experiments is about 1.5x1 0 4 N/m. 


Experiments For Flexible Rotor 

Critical speed map . - Figure 9 shows the critical speed map of the rotor 
model shown in figure 2. In the case shown, both stations 2 and 6 are sup- 
ported by the equivalent spring K. For zero bearing stiffness, there are two 
rigid modes at 0 Hz, the first bending mode at 20 Hz, the second at 86 Hz, and 
the third at 130 Hz. The computed free-free data have good agreement with the 

measured free-free test data shown in figure 10. From this, it seems that the 

mathematical model is correct. The actual critical speeds are estimated (shown 
by dot-dash line in fig. 9) using the measured spring constant. 

Frequency response of flexible rotor - magnetic bearing system . - Figure 11 
shows the measured frequency responses of the FFT analyzer based on an impulse 
response using the same conditions as in figure 8. Figure 12 shows measured 
mode shapes. The first mode at 10.40 Hz is omitted because it is similar to the 
second mode. The modal damping ratio at the first mode was C, = 15.04 percent. 
The measured critical speeds up to the fifth in figures 11 and 12 coincide with 

the computed critical speeds shown in figure 9. The resonance peaks from the 

first through the third are not clearly recognized because of their high modal 
damping ratios as shown in figure 12. 

Eigenvalue analysis . - In this section, we determine the eigenvalues by 
using equation (24). Table I lists the eigenvalues with collocation in the 
ideal case, and table II shows the eigenvalues without collocation in the case 
of this experiment. Only complex eigenvalues are listed in the tables and the 
real roots are omitted. The parameters and feedback gains used for the eigen- 
value analysis are the same as those used in the previous section. The results 
in table II almost coincide with the results of the critical speeds shown in 
figures 9 and 11. The fourth mode is unstable in table II. The measured modal 
damping ratio at 123.13 Hz (as shown in fig. 12) is too small. Moreover, we 
cannot pass the fourth critical speed (as shown in fig. 13) because the ampli- 
tudes at locations near the magnetic bearings are growing at speeds greater than 
5500 rpm. Therefore, the results of the eigenvalue analysis of equation (24) are 
reasonable. It is possible to design the control system of a flexible rotor - 
magnetic bearing system by using the eigenvalue analysis of equation (24). It 
is clear that every mode is stable with collocation. However, the roots of the 
higher order modes approach the imaginary axis as shown in table I. Table III 
shows the improvement of the higher modes with collocations in addition to new 
velocity feedback loops. The difficulty of passing the higher order critical 
speeds with collocation is a very important result. Namely, it points out that 
the usual design method for a rigid rotor has a natural limitation for control 
of flexible modes. New active control loops are needed for the higher flexible 
modes even with collocation. 

Unbalance response . - Figure 13 shows the unbalance responses for this 
experiment for the conditions indicated previously. The results show that both 
the horizontal and the vertical responses have similar characteristics. There 
are only low amplitude resonances up to 5500 rpm. The rotor should have passed 
the fourth critical speed by 5500 rpm since the fourth critical speed exists at 
84 Hz, 5040 rpm as shown in figure 12. Some results indicate that this is 
caused by gyroscopic effects. 
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Active vibration control . - Improving the damping characteristics at the 

fourth mode is needed in addition to considering the velocity feedbacks f^ 

and f^ . Namely, the vibration velocities measured at locations 3 and 5 are 

used as the control inputs at locations 2 and 6. The results of the eigenvalue 
analysis for this are shown in figure 14. The measured frequency responses by 
an FFT analyzer are shown in figures 11(b) and (d). The analysis and experi- 
ments show that the characteristics of the fourth mode are improved by active 
velocity feedback loops. It is possible to pass the fourth critical speed if 
the dynamic characteristics are as shown in figures 11(b) and (d). However, it 
is impossible to pass the fourth critical speed in the actual test rig because 
of gyroscopic effects. We must design the active vibration control system tak- 
ing into account the gyroscopic effects of matrix [G] in equation (22). The 
higher order critical speeds can be passed if the flexible rotor - magnetic 
bearing control system is designed using the stable eigenvalues of equation (24) 
that take into account gyroscopic effects. 


CONCLUSIONS 

This paper deals with vibration problems and control methods without gyro- 
scopic effects for a flexible rotor supported by radial magnetic bearings with 
a four-degree-of-freedom system. The conclusions obtained are as follows: 

(1) The eigenvalue problem, taking into account the interaction between 
the mechanical system of a flexible rotor and the electrical system of magnetic 
bearings and the controller, can be formulated as shown herein. The control 
system of a flexible rotor - magnetic bearing system was designed by an eigen- 
value analysis. 

(2) If the control system is designed for a rigid rotor, the flexible 
modes do not become unstable with collocation. 

(3) New active control loops must be added for the higher flexible modes 
even with collocation because the damping characteristics degenerate for the 
higher order modes. It is possible to stabilize for the case of no collocation 
by means of this method. 

(4) It is not possible to pass the critical speed in the actual test rig 
even if frequency responses by a FFT analyzer are satisfactory mainly because 
of gyroscopic effects with rotation. 
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TABLE I. - EIGENVALUES WITH 
COLLOCATION 


Real 

Imaginary 

Frequency, 

Hz 

-0.22273D-2 

0.361 28D+5 

5750.0 

— . 22359D— 1 

. 1 5725D+5 

2502.7 

- . 70242D-1 

. 12784D+5 

2034.6 

-.5751 9D+0 

. 59381 D+4 

945.1 

-.89774D+1 

.291 59D+4 

464.1 

-.22630D+2 

.20626D+4 

328.3 

-.64734D+2 

.94692D+3 

150.7 

-.47460D+2 

.66731D+3 

106.2 

-.82129D+2 

. 1 5805D+3 

25.2 

-.94146D+2 

. 1 5287D+3 

24.3 

-.49705D+1 

. 12235D+3 

19.5 


TABLE II. - EIGENVALUES WITHOUT 
COLLOCATION 


Real 

Imagi nary 

Frequency, 

Hz 

0.16181 D+2 

0.201 01 D+4 

319.9 

.74842D+2 

.86925D+3 

138.3 

.71 604D+2 

.62756D+3 

99.9 

-.111 1 1D+3 

. 1 7393D+3 

27.7 

-. 12722D+3 

. 1 5032D+3 

23.9 

-.65155D+1 

. 10677D+3 

17.0 


TABLE III. - EIGENVALUES WITH COLLOCATION 
IN ADDITION TO VELOCITY FEEDBACK 


Real 

Imaginary 

Frequency, 

Hz 

-0.4691 1D+3 

0.11 555D+4 

183.9 

-.59602D+3 

.91 103D+3 

145.0 

-.1 4625D+3 

. 36865D+3 

58.7 

-.28629D+3 

. 3451 5D+3 

54.9 

-.46764D+0 

. 1 3729D+3 

21.9 
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FIG. 5.- THREE-DIMENSIONAL STABILITY CRITERION CHART IN 
HORIZONTAL DIRECTION. 








FIG. 10. - MEASURED FREE-FREE TEST DATA. 
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(C) DISK 5. 

FIG. 13. - UNBALANCE RESPONSES. 
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